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Abstract 

We study the existence theory for parabolic variational inequalities in weighted 
L 2 spaces with respect to excessive measures associated with a transition semigroup. 
We characterize the value function of optimal stopping problems for finite and infinite 
dimensional diffusions as a generalized solution of such a variational inequality. The 
weighted L 2 setting allows us to cover some singular cases, such as optimal stopping 
for stochastic equations with degenerate diffusion coefficient. As an application of the 
theory, we consider the pricing of American-style contingent claims. Among others, we 
treat the cases of assets with stochastic volatility and with path-dependent payoffs. 

Keywords: Variational inequalities, excessive measures, Kolmogorov operators, m- 
accretive operators. 
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1 Introduction 

The aim of this work is to study a general class of parabolic variational inequalities in 
Hilbert spaces with suitably chosen reference measures. In particular, our motivation 
comes from the connection between American option pricing in mathematical finance and 
variational inequalities. It is well known by the classical works of Bensoussan [5j and 
Karatzas |18| that the price of an American contingent claim is the solution of an optimal 
stopping problem, whose value function can be determined, in many cases, solving an asso- 
ciated variational inequality (see e.g. [15] for the classical theory and [T7] for connections 
with American options). 

In this paper we study variational inequalities associated to finite and infinite dimen- 
sional diffusion processes in I? spaces with respect to suitably chosen measures. In par- 
ticular, denoting by L the Kolmogorov operator associated to a diffusion AT on a Hilbert 
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space H, we shall choose a probability measure \i that is (infinitesimally) excessive for L, 
i.e. that satisfies L*/x < lo/j, for some u£l (see below for precise statements). An appro- 
priate choice of reference measure is essential in the infinite dimensional case, as there is no 
analog of the Lebesgue measure, and turns out to be useful also in the finite dimensional 
case to overcome certain limitations of the classical theory. In particular, we can relax the 
usual nondegeneracy assumptions on the diffusion coefficient (or on the volatility, using the 
language of mathematical finance), which is usually assumed in the "traditional" approach 
of studying variational inequalities in Sobolev spaces w.r.t. Lebesgue measure (see [6], 
|17j). This allows us, for instance, to characterize the price of American contingent claims 
on assets with degenerate or stochastic volatility as the solution of a variational inequality. 
Similarly, we can treat path-dependent derivatives, as well as claims on assets with certain 
non-Mar kovian price evolutions, using the infinite dimensional theory. We would like to 
mention that Zabczyk [25] already considered variational inequalities (called there Bellman 
inclusions) in weighted spaces with respect to excessive measures, including specific for- 
mulas for excessive measures and applications to American option pricing. However, some 
of our results on existence of solutions for the associated variatonal inequalities are more 
general (our assumptions on the payoff function are weaker, we allow time-dependent pay- 
offs), and we explicitly construct a reference excessive measure in many cases of interest. 
Let us also recall that a study of diffusion operators in IP spaces with respect to invariant 
measures (i.e. measures \x such that L*fi = 0) has been initiated in |23| . 

The main tool we rely on to study the above mentioned optimal stopping problems is 
the general theory of maximal monotone operators in Hilbert spaces. However, we need 
some extensions of the classical results, which are developed below and seem to be new. 
In particular, we establish abstract existence results for variational inequalities associated 
to the Kolmogorov operator of finite and infinite dimensional diffusions (on these lines see 
also [4] and [2l]). 

Variational inequalities connected to optimal stopping problems in finance have also 
been studied in the framework of viscosity solutions, see e.g. [22], [14] . In particular in 
the latter paper the authors consider the problem of optimal stopping in Hilbert space 
and as an application they price American interest rate contingent claims in the Goldys- 
Musiela-Sondermann model. Using the approach of maximal monotone operators, at the 
expense of imposing only very mild additional assumptions on the payoff functions, we are 
able to obtain more regular solutions, which also have the attractive feature of being the 
limit of iterative schemes that can be implemented numerically. Moreover, the additional 
conditions on the payoff function we need are satisfied in essentially all situations of interest 
in option pricing. 

The paper is organized as follows: in section [2] we prove two general existence results for 
the obstacle problem in Hilbert spaces. In section[3]we relate these results with the optimal 
stopping problem in Hilbert space. Applications to the pricing of American contingent 
claims are given in section 4. 
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2 Abstract existence results 



Let us first introduce some notation and definitions. Given any Hilbert space E, we shall 
always denote by | • \e its norm and by (•, ■) E its scalar product. Moreover, we define 
C([0,T], E) as the space of .E-valued continuous functions on [0, T], and W 1,p ([0,T], E), 
1 < p < do, as the space of absolutely continuous functions (p : [0, T] — > E with G 
L p ([0, T], -E). The space of Schwarz' distributions on a domain S C will be denote by 
Similarly, W ,S ' P (S) stands for the set of functions 4> : S — > 1R that are in L P (H) 
together with their (distributional) derivatives of order up to s. Finally, (j) G W£?(S) if 
4>C £ Ty s,p for all C £ Cc°(^)i the space of infinitely differ entiable functions on H with 
compact support. 

Let H be a Hilbert space and [i be a probability measure on i7. Denote by W the Hilbert 
space L 2 (H,fi). Let (-P^)t>o be a strongly continuous semigroup on 7i with infinitesimal 
generator — N. We shall assume that 

\Ptt\n < e wt \<j>\ H Vt > 0, g H, 
where u; G K. Then iV is cu-m-accretive in i.e. 

<iV<M>« > V<peD(N) 

and i?(AJ + iV) = for all A > u, where D(-) and R(-) denote domain and range, 
respectively. Let g G 7i be a given function and define the closed convex subset of TL 

Kg = {(j) £H : > g /x-a.e.}. 

The normal cone to tC g at is defined by 

M g {cj)) = {z G H : J z(<f>- ip)dfjL > W G /C}, 

or equivalently 

= |z G W : = if > g(x), z(x) < if <£(:c) = g(x), /i-a.e. j. 
We are going to study the parabolic variational inequality 

V(0) = <A), 

where ^ G H and / G L 2 ([0,T],ft) are given. 

By a strong solution of (|2.ip we mean an absolutely continuous function (/? : [0, T] 7i 
which satisfies (12. ip a.e. on (0, T). A function G C([0, T],H) is said to be a generalized 
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solution of (|2.ip if there exist sequences {</?o} C Tt, {f n } C L 2 ([0, T], W) and {(/?«.} C 
C([0,T],TC) such that, for all n, ip n is a strong solution of 

^(t)+N<p(t)+Mg(<p(t))Bf n (t) 

a.e. on (0, T) with initial condition 99(0) = <^g, an d fn — ► ¥> in C([0,T],7Y) as n — > 00. 

In order to establish existence of a solution for equation (|2.ip we are going to apply 
the general theory of existence for Cauchy problems in Hilbert spaces associated with 
nonlinear maximal monotone operators (see e.g. [2], [3], [TO]). We recall that the nonlinear 
(multivalued) operator A : D{A) C 7i — > 7i is said to be maximal monotone (or equivalently 
m-accretive) if (y% — 1/2, %i — ^2) > for all yi £ Axi, i = 1,2, and i?(I + vl) = H. The 
operator ^4 is said to be w-m-accretive if XI + A is m-accretive for all A > u. If A is 
cj-m-accretive we set (Yosida approximation) 

1 1 

yl A n = -(u - (I + XA^u), ueH, 0<A<-. 
A w 

Recall that A\ is Lipschitz and -accretive on "H, i.e. 

(A x u - A x v,u- v) H > - - ^ Xuj W - v\h- 

Moreover, recalling that N is w-m-accretive, we have the following result. 

Theorem 2.1. Assume that Pf is positivity preserving (that is Ptip > for all if > 
fi-a.e.) and 

\(N x g) + \n<C VA€(0,l/w). (2.2) 
T/ien the operator N + M g with the domain D{N) Pi /C s is u -m-accretive in H. 

Proof. It is easily seen that iV + Af g + ul is accretive. In order to prove m-accretivity, let 
us fix / £ 7i and consider the equation 

a f\ + N x <p\ + Af g ((p\) B f, (2.3) 

which admits a unique solution for a > u>/(l — Acj), because the operator N\ + M g + al is 
m-accretive for a > oj/(1 — Xu). We are going to show that, as A — > 0, <p\ — * tp strongly in 
H to a solution ip of 

aip + N<p + N g (<p) 3 f. (2.4) 

Let us rewrite ()2.3[) as 

a^A + N x tpx +N K {^x) Bf-ag- N x g, (2.5) 

where tp\ = <p\ — g, K, = {ip £ H : ^ > /z-a.e.}, and jVjc is the normal cone to /C. Setting 
7]\ £ N/c(i/j\) and multiplying both sides of (|2.5p by r/> we have 

a (V>A, »?a)« + + (NaV'a + N x g, rjx) n = (f - ag, r]\} n . (2.6) 
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Since {ip\,Vx)n — (by definition of J\ffc) and {Nxipx^x)^ > (in fact (/ + XN) l K. C K. 
because Pt is positivity preserving), (|2.6f) yields 

\v\\n + (N\g,i]\) H < (f-ag,rj\) H . (2.7) 

On the other hand, we have {N\g,r]\) n > ((N\g) + , because r]\ G Nic{ip\) implies 

that {rj\, (j))^ < if 4> > /i-a.e.. Then by (|2.7h and assumption (12.21) we obtain 

|t/a|w < 1/ " ag\ H + |(AT A5 r)+| w < C VA G (O,^ 1 ). 

Moreover, ()2.5[) implies that 

I^aIw < 1/ - otg\ n VAg^,^- 1 ). 

Therefore = ipx + 5} and {r/^} are bounded in Ti., and so is {N\ip\}. This implies by 
standard arguments that {fx} is Cauchy in 7i, so we have that on a subsequence, again 
denoted by A, 

ipx — ► 93 strongly in 
N\(<P\) -> £ weakly in H, 
?7a — ► V weakly in H, 

as A — > 0. Since t?a G -Xfg(fx) and A/" 9 is maximal monotone, we have r] G N g (ip) and, 
similarly, £ = iV<^. Hence y> is a solution of (|2.4j) . as required. □ 

Remark 2.2. If Pt is the transition semigroup associated to a Markov stochastic process 
X, then is automatically positivity preserving. Assumption (|2.2p holds in particular if 
g G D(iV) or (/ + XN)- 1 g > g for all A G (0, 1/u). 

Remark 2.3. Denoting by N* the dual of N, the operator N has a natural extension 
from H to (D(N*))' defined by JV«(yi) = it(iVV) for all ip G D(JV*) and u G W. Then 
as A — > one has iVxp — > Ng weakly in (D(N*))' and if it happens that Ng belongs to a 
lattice subspace, then condition (|2.2[) simply means that (Ng) + G "H. This is the case in 
spaces L 2 (E), H C M n , where usually A^g is a measure on H (see e.g. [9]). 

Remark 2.4. Theorem 12.11 remains true if we replace assumption (|2.2p by 

\\{g-Ptg) + \n<C Vie (0,1). (2.8) 

The proof follows along completely similar lines. 

By the general theory of Cauchy problems associated with nonlinear m-accretive oper- 
ators (see e.g. [2], [3], [TO]) we obtain the following result. 
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Theorem 2.5. Assume that the hypotheses of Theorem \2.1\ are satisfied. Let (fo G D(N)D 
Kg and f G Ty 1,:L ([0, T];TL). Then there exists a unique strong solution tp G VF 1,oo ([0, T];H)n 
L°°([0, T]; D(N)) of the Cauchy problem h2. 1\) . Moreover the function t \— > ip(t) is right- 
differ entiahle and 

^<p(t) + ct>(t) = 0, te[0,T), 

where 

(Ntp(t)-f(t) /i-a-e. in {<p(t,x) > g(x)} 

4>(t) = s 

{ (Ncp(t) - f(t))+ /i-a.e. in {<p(t,x) = g(x)}. 

If (fo G Kg and f G L 2 ([0,T];TL) then equation \2. 1)) has a unique generalized solution 
p G C([0,T],H), <p(t) G Kg for almost all t G [0,T]. 

We shall see later (see Theorem 12.81 below) that the generalized solution satisfies (|2.ip 
in a more precise sense. 

Remark 2.6. By the general theory of Cauchy problems for nonlinear accretive oper- 
ators (see [2], [3], [10]) one knows that the solution <p{t) given by Theorem 12.51 can be 
approximated as h —* by the solution {pi}^i of the finite difference scheme 

tp l+ i + hNip i+1 + hM g ((p i+ i) 3 fi + tpi, i = 0, 1, . . . , N h , 

where hNh = T and fo = j}l h f{t) dt. Equivalently, 

f <pi +1 = (I + hNyHfo + w), if (/ + hNy^fo + > g, 
\ (fi> g Vi 

2.1 Time-dependent obstacle 

We shall consider the case where the obstacle function g depends also on time. In particular, 
we shall assume that 

gGW l ^([0,T],7i) (2.9) 

f \{N x g)+\ 2 H dt<C VAG^u,- 1 ). (2.10) 
J o 

Let gx = (I + \N)~ 1 g and consider the approximating equation 

^(t)+N(px(t)+gx(t)-g(t)) +M g{t) {px{t)) b f(t) (2.11) 

on (0, T) with initial condition (p(0) = p>o, and ipo > g(0), f G L 2 ([0,T],Tt). Equivalently, 
setting ipx = <px ~ 9, we get 

f + +M K (Mt)) B f(t) - f t (t)-N x g(t), 

MO) = <p o -g(0) G/C. 
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In order to work with strong solutions of equation (|2.12p . we shall assume, without any 
loss of generality, that / G W /1 ' 1 ([0,T],H), % G W 1 ' 1 ([0,T],?i), and p -s(0) G ICPiD(N). 
This can be achieved in the argument which follows by taking smooth approximations of 
/, g and ipo . Then equation (12.12R has a unique strong solution tp x G W 1,oo ([0,T],H) n 
L°°([0, T], D(N)) by standard existence results for Cauchy problems because, as seen ear- 
lier, N + Mjc is w-m-accretive. Moreover, multiplying both sides of (|2.12p by r]\(t) G 
J\fjc(i/j\(t)) and taking into account that Pt is positivity preserving and 



(N^ x , Vx ) n dt>0, £ (^t±( t ), Vx (t)^ dt = VXg^lu- 1 ) 



rT 

io " Jo \ ai In 

arguing as in the proof of Theorem l2.lt we get the following a priori estimates: 

IvaWIm < c ViG[o,r], (2.13) 

T 

\Vx{t)\ndt < C, (2.14) 

o 

for all A G (0, uj^ 1 ). Hence on a subsequence, again denoted by A, we have 

cp x — ► (p weakly* in L°°([Q,T\,H) 
r]x -> rj weakly in L 2 ([0,T},H) 

as A — > 0. Moreover, ip : [0, T] — ► 7i is weakly continuous and 

^(t)+^(i) + # = /(t) (2.15) 
almost everywhere in [0, T] with initial condition <p(0) = <po in mild sense, i.e., 

<p(jt) + f e- N{ - t - s \{s) ds = e- m <p + f e~ N ^f(s) ds 
Jo Jo 

for almost all t G [0,T]. The latter follows by letting A — > into the equation 



VA(t)+PA(*)-p(t)+ / e~ N ^( Vx (s)-f(s)-g' x (s)+g'(s))ds 

= e- Nt (p +g x (0)-g(0)). 



(2.16) 



Taking into account that, as A —* 0, <?a(0 ~~ ► <?(i) strongly in W on [0, T] and g' x — g' = 
(I + XN)^ 1 g' — g' — ► strongly in L 2 ([0, T], W), we obtain the desired equation. In particular 
it follows that (f X (t) — > weakly in 7i for i G [0,T]. We are going to show that 
rj(t) G N g (<p(t)) a.e. on [0,T]. To this purpose it suffices to show that 

rp rp 

limsup / e*(ri x (t),<px(t)) H dt< f (r?(i) , <p(t)) n dt (2.17) 

A^0 JO JO 
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for some real number 7. We shall prove that (|2.17p holds with 7 = — 2u. To this end we 
set N u = N + ujI (note that N u is m-accretive in H) and, rewriting equation (|2.1ip as 

^(^a + g\ - g) + N u (<p x + g\ - g) + v\ - + gx - g) = f + g'x - g' , 

we may equivalently write (|2.16p as 

e-^bx® + 9x(t) ~ g(t)) + f e- N ^e-^( Vx (s) - f(s) - g' x (s) + g'(s)) ds 

Jo 

= e- N « t (Lp +gx(0) -g(0)) VtG(0,T). 

This yields 

f 

Jo 

- jf (e-^rjxit), J* e-W-'h-"'^)} ds 
+ f (vx(t),e- N ^(ipo + gx(0) - g(0)) - e^\g x {t) - g(t))) H dt 



-2ut 



{vx(t),fx(t))n dt 







r / e""* 



H 

(rjx(t),J* e- N ^e-^(f(s) + g' x (s) - g'(s)) ds} dt 



10 \ jo in 
Then letting A —* we obtain 



A-»0 JO 



limsup / e- 2ujt (nx(t), <p\(t)) n dt < 



liniinf^ (e-^ri^t), J e^^s) ds ) dt (2.18) 

T 

cut 



+ J e- ut (^( y t),e- N - t i P o + e~ N ^ t ' s) e- UJS f(s)ds^ dt. 

Consider the function 

F(y) = £ (v(t)>J* e- N " {t - s) y(s) da} dt, y G L 2 ([0, T],H), 

which is continuous and convex on L 2 ([0,T],H.) (the latter is an easy consequence of the 
fact that N u is accretive). Hence F is weakly lower semicontinuous and therefore 

liminf F(e~^m) > F(e -£J V). 

A^O 
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Substituting this expression into (|2,18p we find that 

- 2 - t {r lx {t)^ x {t)) H dt< 



lim sup 

A^O JO 



+ 



r 



T 



—uit 



dt 



— Ujt 



-Nut, 



n 

-Nu>(t—s)—ujs 



e- UJS f{s)ds) dt 



(2.19) 



n 



-2u>t 



77(f), e JV "Vo + 
(v(t),<p(t)) H dt. 

The latter follows by equation dtp/dt + Nip + r/ = f, or equivalently 

- wt f(t). 

Hence T)(t) G N g ((p(t)) for all t G (0,T) as claimed. 

Definition 2.7. A function ip G C([0, T],H) is said to 6e a miW solution of 



%{e^<p(t))+N u (e-^<p(t))+ri(t)e-" t 



^(t)+Nip(t)+M g{t) ( V (t))3f(t) 



(2.20) 



on [0, T] with initial condition (p(0) = <po if v{t) > g(t) \i-a.e. for almost all t G [0, T] and 
there exists rj G L 2 ([0,T],Tt) with rj(t) G N g M(ip(t)) for almost all t G [0,T], such that 



-N(t-s) 



rj(s) ds 



-JV(t-s) 



7(s)ds 



(2.21) 



/or aHt G [0,T]. 

Theorem 2.8. Assume that Pt is positivity preserving and &2.9\) . 112. 1 U\) hold. Let </?q G 
¥>o > 5(0) and f G L 2 ([0, T], W). Then \2. 20\) has a unique mild solution. Moreover, the 
map (<po,f) ^ ip is Lip schitz from H x L 2 ([0,T],Tt) to C([0,T],H). 

Proof. Existence was proved above. Uniqueness as well as as continuous dependence on 
data follows by (I2.2ip taking into account that r](t) G N g (t){<p{t)) for almost all t G [0, T] 
and 



(r)(t), j\- N{t ~ s) V (s)ds^ dt > 



-N(t-s) 



r](s) ds 



dt 



H 



T 2 
p -N{t-s). 



'rj{s) ds 



H 



□ 



It is worth emphasizing that in the case where g is time-dependent the "mild" solution 
provided by Theorem 12.81 is a generalized solution in the sense of Theorem 12.51 However, 
even in this case Theorem 12.81 is not directly implied by Theorem 12.51 
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3 Variational inequalities and optimal stopping problems 

Let H be a real Hilbert space with inner product (•,•) and norm | • |, and (Q,,J-,F = 
(J~t)t>o,W) a filtered probability space satisfying the usual conditions, on which an H- 
valued Wiener process (adapted to F) with covariance operator Q is defined. Let X be the 
process generated by the stochastic differential equation 

dX(s) = b(X(s)) ds + a(X(s)) dW(s) (3.1) 

on s £ [t, T] with initial condition X{t) = x, where b : H — > H and a : H — > L(H, H) are 
such that (|3.ip admits a unique solution that is strong Markov. Define the value function 
v(t, x) of an optimal stopping problem for X as 

v(t, x) = sup E t>x \e-^g(r, X(r)) + f e~^f(s, X(s)) ds] , (3.2) 
rem 1 Jt J 

where 9JI is the family of all F-stopping times such that r £ [t,T] F-a.s., and 

iP(t, s) = J c{X r ) dr \/t<s<T, 

where c : H — > R + is a given discount function (which we also assume to be bounded, for 
simplicity). Exact conditions on g and / will be specified below. The function v is formally 
the solution of the backward variational inequality 

du 
dt 

in (0, T) x H with terminal condition u(T,x) = g(T,x), where 

L O = \ ^[{cyQ^icjQ^yD 2 ^ + (b(x),D<l>) , G £>(L ) = C 2 b {H). (3.4) 

More precisely, denoting by fi an excessive measure of the transition semigroup P% generated 
by the process X, we have that v is the solution of the variational inequality 

du 

— +Lu-cu-jV g(t) (u)Bf (3.5) 

in (0, T) with terminal condition u(T) = g(T), where L is the infinitesimal generator of 
P t . In many situations of interest L = Lq, the closure of Lq in L 2 (H,^l). Before giving a 
simple sufficient condition for this to hold, let us define precisely excessive measures. 

Definition 3.1. Let Pt be a strongly continuous semigroup on L 2 (H,fi), where fi is a 
probability measure on H . The measure [i is called excessive for Pt if there exists uo > 
such that 

[ Ptf dfi < e w * / fdfi Vt > 
for all bounded Borel functions f with f > fi-a.e.. 



+ L u - cu- N g{t ) (u) B f (3.3) 
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We have then the following result. 

Lemma 3.2. Let Pt the semigroup generated by X, and let [i be an excessive measure for 
P t on H. Moreover, let b G C 2 {H) n L 2 {H, y), a G C 2 (H, L(H, H)), and 

\Db(x)\ H + \Da(x)\ L(H<H) < C (3.6) 



for all x G H. Then —Lq is uj-accretive and L is the closure in L (H,fi) of Lq defined on 
D(L ) = C 2 (H). 

Proof. The argument is similar to that used in [12] for similar problems, so it will be 
sketched only. Fix h G C 2 (H) and consider the equation (XI — Lq)<p = h, or equivalently 



poo 

ip(x)=E 0tX e- xt h(X(t))dt, X>lo. (3.7) 
Jo 

It is readily seen that ip G C 2 (H) and, by Ito's formula, (A — Lo)ip = h in H. Since — Lq 
is closable and w-accretive, and R(X — Lq) is dense in L 2 (H,/j 1 ), we infer that Lq coincides 
with L. □ 

Note also that since the measure \i is w-excessive for Pt we have J H Lf dy < to f H fdy, 
which implies that L is o;-dissipative in L 2 (H, fj,). In the sequel, for convenience of notation, 
we shall set N = — L + cl. 

We shall further assume that g(t, x) is continuously differentiable with respect to t, 
Lipschitz in x, and 

sup / (\D t g(t,x)\ 2 + \D x g(t,x)\ 2 )n(dx) < oo, (3.8) 
te(o,T) JH 

Tx[{aQ l ' 2 ){aQ l ' 2 YD 2 xx g] > on (0,T) x H. (3.9) 



If H is a finite dimensional space, the inequality (|3.9p must be interpreted in the sense of 
distributions (i.e. of measures). In the general situation treated here the exact meaning of 
(|3.9p is the following: there exists a sequence {g e (t)} C C 2 (H) such that 

sup [ (\D t9e (t,x)\ 2 + \D x9e (t,x)\ 2 )fi(dx) < C Ve>0, 
te(o,r) JH 

TrKaQ^Xo-Q 1 / 2 )* D 2 xx g £ {t,x)} > Ve > 0, t > 0, x G H, 

9e{t)^g{t) in L 2 (H,y) Vi > 0. 

It turns out that under assumption f|3.9|) g satisfies condition (|2.10p . Here is the argument: 
for each A > we have (N\g) + = lim £ ^Q(N\g £ ) + in L 2 (H,y). On the other hand, N\g e = 
N(I + XN)- 1 g e and by ^M) we see that 



Tr 



(aQ^iaQ^YDl^I + XN)-^,] 



> on H 
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because (I + XN) 1 leaves invariant the cone of nonnegative functions (by the positivity 
preserving property of Pj). Hence 



(N x g £ 



< 



(6, + AJV)" 1 ^ 



L 2 (ff,/x) 



<C VAs(0,^ 1 ),e>0 



because 6 G L 2 (H,fi). This implies (|2.10p as claimed. 



Proposition 3.3. Assume that f G L 2 ([0, T],L 2 (H, fi) nC([0,T],C b (H)) and that condi- 
tions h2.9\) . \3. 6\) and \3. 9\) hold. Furthermore, assume that the law of X(s) is absolutely 
continuous with respect to /j for all s G [t,T]. Then there exists a unique mild solution 
u G C([0, T]; L 2 (H, fj,)) of the variational inequality h3. 5\) . Moreover, u coincides fi-a.e. 
with the value function v defined in &3.2\) . 



Proof. Existence and uniqueness for (|3.5p follows by Proposition 12.81 In the remaining of 
the proof we shall limit ourselves to the case / = 0. This is done only for simplicity, as 
the reasoning is identical in the more general case / 7^ 0. By definition of mild solution 
there exists r? G L 2 ([0, T],L 2 (H, /x)) such that rj(t) G N g{t) (u(t)) for all t G [0,T] and the 
following equation is satisfied (in mild sense) for all s G (0,T), with terminal condition 
u(T)=g(T): 

fin 

(s)-Nu(s)=ri(s), (3.10) 



i.e., 



u(t, x) 



ds 



Rs-tVi s i x ) ds + R T -t u ( T , x) Vi < t < T, /j-a.e. x G H, 



(3.11) 



where Rt is the transition semigroup generated by — N, or equivalently the following 
Feynman-Kac semigroup associated with the stochastic differential equation (|3.1f) : 



R t cf>(x) = E 0tX \e- J"o ds 0(X(t)) 



GL 2 (H,fi). 



Let us set Ht = [t, T] x H and define the measure \xt = Leb x \x on Pr, where Leb stands 
for one-dimensional Lebesgue measure. Recalling that u(s,x) > g(s,x) for all s G [i, T], 
/j-a.e. x £ H, we can obtain a version of u, still denoted by u, such that ii(s, ac) > g(s, x) 
for all (s,x) G Pt- Recalling that r/(s, •) G L 2 (P, /j) a.e. s G [i,T], equation (|3.1ip yields 



u(t, ac) = Et x 



-^r]{s, X(s)) ds + e-^(*' r ^(r, X(t)) 



(3.12) 



for every stopping time r G [t, T], for all t G [0, T] and /i-a.e. x £ H. In fact, let us 
consider a sequence {r/ £ } C C 1 ([0, T],W) such that n £ — ► 77 in L 2 ([0, T], 7i). Then equation 
dSHOD, with % replacing 77, admits a solution u £ G C 1 ([0, f],W)nC([0, T], P(iV)) such that 
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u £ — » u in C([0,T],W) as e — > 0. Recalling that N = —Lq + cJ, there exists a sequence 
{u>J C C 1 ([0,r|,W) nC([0,T],Cf(iJ)) such that 

\u £ (t) - w £ (t)\ n < e 
\Nu £ (t)-(-L + cI)w £ {t)\ H < e 
. du £ , . dw £ , . . 

I-5-W-^WIk < . 

for all i G [0, T]. Therefore we have 

dw, 



(t) - (Lq - cl)w e (t) = fj £ VtG[0,T], 



(it 

where % — > 77 in L 2 ([0, T], TC). Then we have 

w £ (t, x) = E M [ [ T - e -^fj £ (s, X(s)) ds + e-^w £ (T, X(r))] (3.13) 



for all stopping times r G [t, T]. We shall now show that (assuming, without loss of 
generality, tp = 0) 

E t>x / fj E (s,X(s))ds^E tiX r](s,X(s))ds 
Jt Jt 

for all i G [0,T] and in L 2 (/7, /j) w.r.t. x. In fact, Tonelli's theorem yields, recalling that 
ji is excessive for Pt, 

' Et >x [ \fj £ (s, X(s)) - V {s, X(s))\ 2 ds »{dx) 



T 

2 



H 

Jo Jh 

rT r 

I P s \f) £ (s,x) - r](s,x)\ 2 /j,(dx)ds 
/o Jh 

— e ^ T I I {ijeisjx) — r](s,x)\ 2 /i(dx)ds —>■ 
Jo Jh 

as e — > 0, because fj £ — > r\ in L 2 ([0,T], L 2 (H, fi)). An analogous argument shows that 
E m w £ (t,X(t)) -» E^-u(t,X(t)) for all i G [0,7] and in L 2 (H,^) w.r.t. z. Therefore, 
passing to a subsequence of e if necessary, we have that the left-hand and right-hand side of 
(|3.13p converge to the left-hand and right-hand side, respectively, of (|3.12[) for all t G [0, T] 
and /x-a.e. x G H. Recalling that 



77(5,2) 



= if u(s,x) > g(s,x) for each s and /j-a.e. x G H, 
< if u(s,x) = g(s,x) for each s and /j-a.e. x £ H, 
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let us define the set 

A = {{s,x) £ H T : rj(s,x) > 0}, 

for which we have (J>t(A) = 0. Using this fact together with the assumption that the law of 
X(s) is absolutely continuous w.r.t. fi for all s £ [t,T], hence that Pt iX ((s, X(s)) £ A) = 0, 
we get 

E tjX j* -e-^ V (s, X(s)) ds > 0. 

Therefore equation (|3.12|) implies that u(t, x) > M tjX [e~^^' T ^ g(r, X(r))] for all stopping 
times r G S9T, hence u(t,x) > v(t,x), for all i G [0,T] and fi-a.e. x £ H. Let us now prove 
that there exists a stopping time f G [i, T] such that u(t,x) = ¥ lttX [e~^^ t,T ^g(f,X(f))], 
which will yield u(t,x) = v(t,x), for all t G [0, T], /i-a.e. x G H. Define the set 

B = {(s, x) G H T : g(s, x) = u(s, x)} 

and the random time 

D B = inf{s >t:(s, X(s)) £ B} AT. 

Since B is a Borel subset of Ht and the process (s,X(s)) is progressive (because it is 
adapted and continuous), the debut theorem (see T.IV.50 in [13]) implies that Db is a 
stopping time. Recalling that u(s,x) > g(s,x) for all s £ [t, Db), we have, reasoning as 
before, rj(s,X(s)) = a.s. for each s £ [t, Db)- Thus, taking f = D B , (|3.12p yields 

u(t, x) = E fiI [e"* ( ' lf) g(T, X(f))\ Vt G [0, T], fi-a.e. x £ H. 

We have thus proved that there exists a version of u such that u(t,x) = v(t,x) for all 
t £ [t,T], fi-a.e. x £ H. The definition of mild solution then implies that u(t,x) = v(t,x) 
for all t £ [t, T] and ^t-a.e. x £ H. □ 

Remark 3.4. The absolute continuity assumption in proposition 13.31 can be difficult to 
verify in general. However, it holds in many cases of interest. In particular, it is auto- 
matically satisfied if the semigroup Pt is irreducible and [i is invariant with respect to Ft- 
Moreover, in the finite dimensional case, if the excessive measure \x is absolutely continu- 
ous with respect to Lebesgue measure and the coefficients of (|3.1j) satisfy an hypoellipticity 
condition, the assumptions of the above proposition are also satisfied. We shall see in the 
next section that /i has full support in all examples considered. Moreover, in the finite di- 
mensional cases, \x can be chosen absolutely continuous with respect to Lebesgue measure. 
Let us also remark that the continuity of the value function has been proved under very 
mild assumptions by Krylov [20], and by Zabczyk [21] in the infinite dimensional case. 

Remark 3.5. Optimal stopping problems in Hilbert spaces and corresponding variational 
inequalities are studied by Gatarek and Switch p3] in the framework of viscosity solutions. 
Their results are applied to pricing interest-rate American options, for which the natural 
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dynamics is infinite dimensional (e.g. when choosing as state variable the forward curve). 
At the expense of assuming (|3.9|) . that is, roughly speaking, a convexity assumption on the 
payoff function g, we obtained here a more regular solution. We would like to remark that 
g is convex in practically all examples of interest arising in option pricing, some of which 
are investigated in the next section. 



4 Pricing of American options 

Let Q be a risk neutral martingale measure, and assume we have n assets whose price- 
per-share X{t) = {Xi(€))i=\ n evolves according to the following Markovian stochastic 
differential equation: 

dX(t) = rX(t)dt + a(X(t))dW(t), X(0) = x>0, t e [0,T], (4.1) 

where r G R + is the risk-free interest rate, W is a M m -valued Wiener process, and a : W 1 — ► 
L(M m , R n ) is the volatility function. Moreover, we assume that a is such that X(t) G M™ for 
all t G [0, T]. The standard assumption (see e.g. [E]) is that aij(X(t)) = Xi{t)dij{X{t)) 
for some a : R n — > L(R m ,R n ). We do not assume, however, that a nor a satisfies a 
uniform nondegeneracy condition. Note that in this situation the market is incomplete, 
even if m = n, and the choice of the risk neutral measure Q is not unique (|19j). 

It is well known that the problem of pricing an American contingent claim with payoff 
function g : R n — * R is equivalent to the optimal stopping problem 

v(t, x) = sup E t>x [e- rT g(X( T ))], (4.2) 
rem 

where 971 is the set of all F-adapted stopping times r G [t, T] and E stands for expectation 
with respect to the measure Q. Denote by P% the transition semigroup associated with 
flUD, i.e. P t f(x) = E , x f(X(t)), f e C b (M. n ), x € R n , and let L be the corresponding 
Kolmogorov operator. A simple calculation based on Ito's formula yields 

L f(x) = ±Tr[a(x)o-*(x)D 2 f(x)] + (rx,Df(x)) Rn , / G C 2 b (R n ). 

By classical results (see e.g. [20J ) , the value function v(t,x) is expected to satisfy the 
following backward variational inequality 

f max ((<9 t + L )v(t,x) - rv(t,x),g(x) - v(t,xfj = 0, (t,x)eQ T 
\ v(T,x) = g(x), ' x G Wl, 

where Qt = [0, T] x R" . 

The classical theory of variational inequalities in Sobolev spaces with respect to Lebesgue 
measure does not apply, however, mainly because the volatility coefficient is degenerate (see 
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|17j). Nonetheless, one might try to study (|4.3[) in spaces of integrable functions with re- 
spect to a suitably chosen measure. The most natural choice would be an (infinitesimally) 
invariant measure for Lq. However, without non-degeneracy conditions for a and with 
r > 0, one may not expect existence of an invariant measure (see e.g. pQ, [7]). Here we 
shall instead solve (|4.3|) in L 2 (R n ,/z), where p is an (infinitesimally) excessive measure for 
Lq, which is also absolutely continuous with respect to Lebesgue measure. 

The backward variational inequality (|4.3p can be equivalently written as the (abstract) 
variational inequality in L 2 (M n ,/x) 

dtv -Nv- M g {v) 3 0, v(T) = g, (4.4) 

where N = —L + rl, with L the generator of Pt (which will often turn out to be the closure 
of Lo), and M g is the normal cone to 

Kg = {</> G L 2 (R n , p) : (J) > g ^-a.s}. 

Lemma 4.1. Assume that 

a£C 2 (R n ), \a(x)\<C{l+\x\), \a Xi \ + \a XiXj \ < C. (4.5) 

Then there exists an excessive probability measure /i of Pt of the form 

= 1 + [^(n+l) dX 

with a > 0. 

Proof. Setting p(x) = 1+ \ x ^ n +i) > we shall check that L^p < cop in W n for some u > 0, 
where Lq is the formal adjoint of Lo, i.e. 

L* p=^Tr[D 2 (aa*p)]-rdiv(xp). 

Assumption (j4.5j) implies, after some computations, that 

sup =: u < oo, 

xei" P 

thus p,(dx) = ap(x) dx, with a~ 1 = f Rn p(x) dx, is a probability measure and satisfies 
LqP < up. This yields 

L Q f dp<u I fdp (4.6) 



for all / G C b 2 (]R n ) with / > 0, and therefore 



Ptfdp<e M fdp 

for all / G C fc 2 (M n ), / > 0. The latter extends by continuity to all / G C b (R n ), f > 0, and 
by density to all bounded Borel / with f > p-a.e.. □ 
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The operator A is w-dissipative in L 2 (M n ,/x). More precisely, we have 

/ (L f)fd^<-\f \(aa*) 1/2 Df\ 2 dy + u [ f 2 dfx V/ G Cf (R n ), 

JR n £ JR n JR n 

as follows by gSD and L (f 2 ) = 2(A)/)/ + ((atr*) 1 ^/^ 
Note also that for each /i G C 2 (M n ) the function 

/>oo 

<p(x)=E 0>x e~ xt h(X(t))dt (4.7) 

JO 

is in C 2 (IR n ) and satisfies the equation 

Xip — Lq^p = h 

in W 1 . Hence R(XI — Lq) is dense in L 2 (M. n , y) and since A is closable, its closure L := Lq 
is w-m-dissipative, i.e. —ujI + L is m-dissipative. Since, by (|4.7p . (XI — L)^ 1 is the resolvent 
of the infinitesimal generator of Pj, we also infer that L is just the infinitesimal generator 
of Pt- We have thus proved the following result. 

Lemma 4.2. The infinitesimal generator of Pt in L 2 (W n ,fi) is L. Moreover one has 
[ (Lf)fdn<-\( \a*Df\ 2 dfi + u; [ f 2 dy 

JR n * JR n JR n 

for all f G L 2 (R n ,y). 

Taking into account that L is the closure (i.e. Friedrichs' extension) of Lq in L 2 (W,y), 
it follows that for each / G D(L) we have 

Lf=^Tr[aa*D 2 f] + (rx,Df) Rn 

in T>'(M. n ), where Df, D 2 f are taken in the sense of distributions. In particular, it follows 
by the previous lemma that 

K) 1/J /6^ 1,2 (R n ,/') ) /e<' c 2 (S) 

for each / G D(L), where S = {x G R n : Tr[aa*](x) > 0}. 

We are now going to apply Theorem 12.11 to the operator N = —L + rl on the set 

K. g = G L 2 (R n ,y) : tp{x) > g{x) /i-a.e.}. 

The function g : W 1 — > M. is assumed to satisfy the following conditions: 

Dg G L°°(IR n ), Tv[aa*D 2 g] G M(R n ), Tx[aa*D 2 g] > in M(R n ), (4.8) 
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where 7W(M n ) is the space of bounded Radon measures on W n . 

Payoff functions that can be covered in this setting include so-called Margrabe options 
(with payoff g(x) = (xj — Xxj) + , for given A > and i ^ j < n) and basket put options. 
We shall focus, as an illustration of the theory, on the latter case, for which 

n n 

g(x 1 ,...,x n ) = (k-^XjxA , ^Aj = l. 

3=1 3=1 

In this case the first two conditions in (14, 8p are obviously satisfied and 

n 

Tr[aa*D 2 g] = [Yl a ii A i A j)^ > 
i,j=l 

in A4(M n ), where a = era* and 5 is the Dirac measure. Moreover, 

n 

Dig(x) = -XiH^r - ^ XjXjJ, i = l,...,n, 

3=1 

where 77 is the Heaviside function, i.e. H(r) = 1 for r > and 77(r) = otherwise. The 
operator N has a natural extension to functions g satisfying the first two conditions in 
(14. 8 p through the formula 

(Ng)(p= j gN*ipdfi V</3 G D(N*). 

In our case one has 

^ n n n 

Ng = ~ 2 ( "U x > x .,) s r Yl X i x i H { k ~Yl X i x i) + rq - 

i,j=l j=l j=l 

Lemma 4.3. Assume that ( f^.5| ) and ft4-8\ ) are verified. Then the operator N + j\f g is 
uj-m- accretive in L 2 (K n ,/x). 

Proof. We only have to prove that \{N\g) + \ L 2^r l ^ is bounded for all A G (0,cj _1 ), as 
required by Theorem 12.11 Set g\ = (7 + XN)~ 1 g, i.e. 

g x + XNg x = g, N x g = Ng x . (4.9) 

Then we have 

A n n 
(1 + Xr)g x (x) - - a ij D ijS\(x) - rJ2 x i D i9\(x) = d(x) 

i,j=X i=l 
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in V(M. n ). As seen earlier, Ng = -\ Tr[aa*D 2 g] - (rx,Dg) + rg in V(W n ) and by as- 
sumption (|4.8j) we have that (Ng) + = (—r(x,Dg) + rg) + (where v + denotes the positive 
part of the measure v). Since Dg £ L°°(W n ,dx) we conclude that (Ng) + £ L 2 (R n ,/x). 
Approximating g by a sequence g e £ D(N) we may assume that 5 £ D(N) and also 
iVc/A £ £>(A0- We set ^ A = iVc/ A and so flM]) yields 

A n n 
(1 + Ar)^A(x) - - a ij(x) D fjMx) -rJ2 x i D iM*) = N 9(x). 



2 

1,3=1 i=l 



Let us set ip\ = ip\ + "0? j with 



A n n 

(l + Ar)^(x)--^a JJJ D^i_ r ^ XiAV; i = (iVff) - 



i,j=l i=l 



, n n 

(1 + Ar)^( x )__ ^ a l3 D}^\ -rJ>A^ = (Ng)~, 

i,j=l i=l 

where the first equation is taken in L 2 (M n , n) and the second in 2?'(R n ). By the maximum 
principle for elliptic equations we infer that ip\ > 0, ip\ > 0, hence tp\ = tp^ and ip\ = . 
This implies that tp't = (Nxg)^ is the solution if)\ of 

4>l + \N4>\ = (Ng) + . 

But the solution of this equation satisfies 

1-1,2 < \( N 9) + \l*(R",h) 

hence {(N\g) + }\ is bounded as claimed. □ 

Applying Corollary 12.51 we obtain the following existence result for the value function 
of the optimal stopping problem, i.e. for the price of the American option. 



Corollary 4.4. Assume that conditions j^.5| ], ( f^.ff| j hold and that g £ D(N) = L 2 (M. n ,fi). 
Then the backward variational inequality associated to the optimal stopping problem J-^.i^ ), 
i.e. 

du 

__ iVu _ A A 9 ( n ) 9 o, a.e. i£(0,T), 
u(T) = g, 

admits a unique generalized (mild) solution u in C([0, T]; L 2 (IR n , fi)). Moreover, if g £ 
D(N), then u £ Vr 1,oo ([0, T]; L 2 (R n , fi)) is the unique strong solution of [JHfy . Further- 
more, if the law of the solution of J^, i| j is absolutely continuous with respect to [i, then the 
value function v coincides with u for all s £ [t,T] and ji-a.e. x £ W 1 . 
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Let us remark that the last assertion of the corollary is included for completeness only, 
as we do not know of any option whose payoff g is smooth enough so that g £ D(N). On 
the other hand, the general case g G D(N) covered in the corollary happens for virtually 
all payoff functions g. Then the solution is just the limit of the following backward finite 
difference scheme: 

Vi = 9 M -i, 0i+i + hN9 l+1 + N g {9 l+1 ) 3 Oi, 9 = v , h = T/M. 

This discretized elliptic variational inequality can be solved via the penalization scheme 

fi? +1 + hN9f +1 - -g)-=9l, i = 0, 1, . . . , M - 1, 

or via the bounded penalization scheme (see e.g. [8]) 

9f +1 + hN9f +1 + gi ^+\~ g =0l + 9i, i = 0, 1, . . . , M - 1, 
£ + ri+l 9\ 

where g± is an arbitrary parameter function. Therefore the characterization of the option 
price given by Corollary 14.41 is also constructive, that is, it is guaranteed to be the unique 
limit of very natural finite difference approximation schemes, that can be implemented 
numerically. A completely analogous remark applies also to the cases treated in the next 
subsections. 

4.1 American options on assets with stochastic volatility 

Consider the following model of asset price dynamics with stochastic volatility under a risk 
neutral measure Q: 

dX{t) = y/V(fjX(t)dWi(t) 

dV(t) = k(9 -V(t))dt + r]^/V{t)dW2{t), 

where W(t) = (Wi(t), Wb(t)) is a 2-dimensional Wiener process with identity covariance 
matrix (the more general case of correlated Wiener processes is completely analogous), k, 
9, 7) are positive constants, and the risk-free interest rate is assumed to be zero. Moreover, 
in order to ensure that V(t) > Q-a.s. for all t G [0, T], we assume that 2k9 > rf (see e.g. 

It is convenient to use the transformation x(t) = log X(t), after which we can write (by 
a simple application of Ito's lemma) 

d x {t) = -V{t)/2dt + y/V(f) dWi(t). 

Define Y(t) = (x(t), V(t)). Then we have 

dY(t) = A(Y(t)) + G(Y(t)) dW(t), (4.11) 
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where A : R 2 B (x,v) i-> (-v/2,k(6 - v)) £ R 2 and G : R 2 3 (x,v) i-> diag(^/v, 77^/i;) G 
L(R 2 , R 2 ). The price of an American contingent claim on X with payoff function g : R — > R 
is the value function v of an optimal stopping problem, namely 

v(t,x,v) = supE L{XtV) [g(Y(T))], (4.12) 
rem 

where g(x,v) = g(e x ) and 9Jt is the set of all stopping times r such that r E [s,T] Q-a.s.. 
The Kolmogorov operator Lq associated to (|4.1ip is given by 

Lof = \vf xx + K, 2 vf vv - l -vf x + k(0 - v)f v , f € CZ(R 2 ), 
and its adjoint Lq takes the form 

LoP = -vp xx + -r] 2 (v P ) vv + -vp x -K((e-v)p) v , P £Ci(R 2 ), (4.13) 

Following the same strategy as above, we look for an excessive measure of the form 

1 



1 + x 2 + v 2 ' 



/i(dx, dv) = ap(x, v) dx dv, p(x, v) 

where o _1 = L 2 p(x, v) dx dv. 

Some calculations involving (|4,13p reveal that 

LZp(x,v) 
sup — — — = OJ < oo, 

i.e. p, is an infinitesimally excessive measure for Lq on 3 = R x R + . Then the transition 
semigroup 

P t f(x,v) = E 0t(x>v) f(x(t), V(t)), f e C 2 (E), 

extends by continuity to L 2 (E, p), and the operator Lq with domain C 2 (E)) is u>-dissipative 
in L 2 (E,p). Arguing as above (see Lemma [372]) . the closure L of Lq is w-m-dissipative in 
L 2 (E, p) and 

JjLf)f dp < -t±l jf v (/2 + /2) dp + U jj 2 dp. 

The operator iV = — L is therefore w-m-accretive and formally one has 

N~g = -± v (e 2x g"{e x )+e x g'{e x )) + l -ve x g'{e x ) 

= -\ve 2x g"{e x ). (4.14) 

The previous expression is of course rigorous if g is smooth and Ng G L 2 (H,//), but in 
general (i.e. for g G L 2 (H,/x)) is has to be interpreted in the sense of distributions on H in 
order to be meaningful. 
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We shall assume that the payoff function g is convex on R, more precisely, 



g" G M(R), g" > 0, (4.15) 

where (R) is the space of finite measures on R. Note that the typical payoff of a put or 
call option is covered by these assumptions. Equation (|4.14p implies that Ng G T>'(M) and 
Ng < in 2?'(R), hence Ng is a negative measure and so the hypotheses of Theorem 12.11 
are met. Thus, defining K, g = {<p G L 2 (E,[i) : cp > g /i-a.e.}, it follow that the operator 
N + M g is cj-m-accretive on H = L 2 (H,/z). This yields 

Corollary 4.5. Assume that j^. 15\ ) holds. Then the backward variational inequality 

du 
~dt 

on Ht = [0, T] x L 2 (H,/i) with terminal condition u(T) = g has a unique generalized (mild) 
solution u G C([0, T], L 2 (H, /x)). Moreover, if g G D(N), then {4-16 ) has a unique strong 



Nu-M g (u)3 (4.16) 



solution u £ W ,oo ([0, T],L (H, fj,)). Furthermore, if the law of the solution of \J^.ll ) 



is 



absolutely continuous with respect to fi, then the value function v defined in QTT^D coincides 
with u for all s G [t,T] and \x-a.e. (x,v) G S. 



4.2 Asian options with American feature 

Let the price process X of a given asset satisfy the following stochastic differential equation, 
under an equivalent martingale measure Q: 

dX = rXdt + a(X)dW(t), X(0) = x. 

Here we consider the problem of pricing a "regularized" Asian options with American 
feature, that is we look for the value function v of the optimal stopping problem 

v(x) = sup E x ( k l — [ X s ds] , (4.17) 

rem V r + o Jo ) 

where k > is the strike price, 5 > is a "small" regularizing term, 9JI is the set of stopping 
times between and T, and E x stands for expectation w.r.t. Q, conditional on X(0) = x. 
The standard Asian payoff corresponds to 5 = 0. Unfortunately we are not able to treat 
with our methods this limiting situation, as it gives rise to a singularity in the obstacle 
function of the associated variational inequality, or, in the approach we shall follow here, 
in the Kolmogorov operator of an associated stochastic system. However, it is clear that 
for small values of 5 the value function v in (14.171) is a good approximation of the option 
price, at least for optimal exercise times that are not of the same order of magnitude of 5. 
Let us define the auxiliary processes 
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and S(t) = t. Then we have 

' dX{t) = rX(t) dt + a(X(t)) dW(t) 

dS(t) = dt 

with initial conditions X(0) = x, Y(0) = 0, s(0) = 0. This system can be equivalently 
written in terms of the vector Z = (X, Y, S) as 

dZ{t) = A{Z{t)) dt + G{Z(t)) dW(t), Z(0) = (x, 0, 0), (4.18) 

where A : R 3 -> R 3 , A : (x,y,s) i-> (rx, (s + - y), 1) and G : R 3 -> L(R,R 3 ) ~ R 3 , 

G(x,y,s) = (cr(x), 0, 0). Therefore (|4.17p is equivalent to 

v(x) = sup E x g(Z(T)), 
rem 

where g : (x, y, s) — > (k — y) + and E x . stands for E^q.o)- 

As in the previous cases, we shall look for an excessive measure of Lq, the Kolmogorov 
operator associated to (|4. 18[) . which is given by 



\a 2 {x)D 2 x J + rxD x f + ^D y f + D.f, f € C 6 2 (M 3 ) 



Then the adjoint of Lq can be formally written as 

L oP = \ D lxie 2 {x)p) ~ rD x (xp) - D v (^-^pj - D sP . 

In analogy to previous cases, some calculations reveal that, under the assumptions (|4.5I) 
on a, there exists a measure p of the type fj,(dx, dy, ds) = p(x, y, s) dx dy ds, 



p(x,y,s) = j— - h2n+ln - j h2n+Ul T a 1= / P( Z ) dz > 

(1 + \x\) zn+1 (l + \x- yiy^+^l + s) z J R 3 

such that LqP < uip for some Arguing as before, we conclude that p is an excessive 

measure for the semigroup Pt generated by the stochastic equation (14.181) . and that L, the 
closure of Lq in I? (R 3 , p) , is the infinitesimal generator of Pt . 

We are now in the setting of section El i.e. we can characterize the option price as the 
(generalized) solution of a suitable variational inequality. Details are left to the reader. 



23 



4.3 Path-dependent American options 

We shall consider a situation where the price dynamics is non-Markovian as it may depend 
on its history, and the payoff function itself is allowed to depend on past prices. We should 
remark, however, that in the present setup we still cannot cover Asian options of the type 
discussed in the previous subsection, with 5 = 0. 

Consider the following price evolution of n assets under a risk-neutral measure Q: 

dX(t)=rX(t) + a(X(t),X s (t))dW(t), 0<t<T 
X(0)=x , X s (0) = x 1 {s), -T<s<0, [ ' y> 

where X s (t) = X(t + s), s £ (— T, 0), W is a standard Wiener process on R n and a : 
R n x L 2 ([— T, 0], R n ) — > L(R n ,R n ). Let us consider an American contingent claim with 
payoff of the type g : R n x L 2 ([-T,0] -> R, whose price is equal to the value function of 
the optimal stopping problem 

v(s,x ,xi) = supE Sj{x0jXl) [e~ rT g{X(T), X s (t))}, 
rem 

where the notation is completely analogous to the previous subsection. An example that 
can be covered by this functional setting is g(xo,x\) = «o5o( 2; o)+02fl , i(^i), with a\, ct2 > 
and go{x ) = (k - x ) + and gi(xi) = {k\ - J° T xi(s) ds) + . 

Let us now rewrite (I4.19P as an infinite dimensional stochastic differential equation on 
the space H = R n x L 2 ([-T, 0], R n ). Define the operator A : D{A) C H -» H as follows: 

A : (x ,xi) i ^ (rx ,x[) 

D(A) = {(z ,xi) eil; x 1 €W 1 ' 2 ((-T,0),R n ), xi(0)=ar o }. 

Setting G(xo,xi) = (a(xo, x±),0), let us consider the stochastic differential equation on H 

dY(t) = AY(t) dt + G(Y(t)) dW(t) (4.20) 

with initial condition Y(Q) = (xq,xi). The evolution equation (|4.20p is equivalent to 
(|4.19j) in the following sense (see [TTj): if X is the unique solution of (|4.19p . then Y(t) = 
(X(t),X s (t) is the solution of (|4.20p . Note that ()4.20p has a unique solution if G is Lipschitz 
on H. Finally, regarding g as a real-valued function defined on H, we are led to study the 
optimal stopping problem in the Hilbert space H 

v(s,x) = sup E S:X [e- rr g(Y(T))}. 
The Kolmogorov operator Lq associated to (|4.20p has the form, on C 2 (H), 

L tf(x ,xi) = -Tr[aa*(x ,x 1 )Dl (p(x ,x 1 )] + (rx ,D X0 (p(x ,xi)) Mn 
+ / {xi(s),D Xl ip(x ,xi(s))) Rn ds. 
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We look for an excessive measure p for Lq of the form p = v\ (g) u 2 , where u\, v 2 are 
probability measures on R n and L 2 {[— T, 0], K n ), respectively. In particular, we choose 



1 



vi(dx ) = p(x )dx , p(x ) = — - — [J-, a=[ — — — ^ efe 

and v 2 a Gaussian measure on L 2 ([— T, 0],M n ). Setting iT = -L 2 ([— T, 0], M n ), a simple 
calculation reveals that 

' L <pdp = - [ dv 2 [ Ti[aa*Dl o (p] dv x + r [ dv 2 [ (x ,D xo ip) Rn dv\ 

H 1 J H if" JH JR" 

+ / dui (xi,D Xl ip) H dv 2 
JR" Jh 

= \ I dv 2 I pD 2 XQ {aa*p)dx -r [ du 2 \ ipD X0 (x p) dx 
1 J H JR n J H il" 

+ / dv x \ (x 1 ,D xl (p) H du 2 . (4.21) 

il™ JHo 



We shall assume that 

a e C 2 (R n x L 2 ([-T,0],lR n )) nLip(M n x L 2 ([-T,0],M n )), 

a{x Q ,x{) < C{\x Q \ + \xx\ho) V(a? , a?i) G (4-22) 

\<rxi{xQ,xx)\ + \a Xi;C:j {xQ,x-i)\<C V(x ,xi) e H, i,j = 1,2. 

Note that these conditions also imply existence and uniqueness of a solution for (|4.19p . 
Taking into account that J H \xi\ 2m dv 2 < 00 and that 

/ {x 1 ,D Xl ip) H dv 2 = - tp(n- {Q~ 1 xi,Xi) H )du 2 
J H () J Ho 

(where Q is the covariance operator of is 2 ), we see by (|4.2ip and (|4.22p that there exists 
u> > such that 

/ L(jipdp<u> / ip dp (4.23) 
JH Jh 

for all Lp E C$(H), p > 0. 

The operator Lq is thus closable and w-dissipative in L 2 (H, p). Moreover, ([4.230 implies 



that 



f (L W )pdp<-\ f \{aa*) l / 2 D Xo p\ 2 dp + uj [ p 2 dp \/ip G C 2 {H). (4.24) 
J h 2 Jh jh 



Since one has, for A > u, 



/•oo 

(\I-L )- 1 <p = E e- xt p(X(t),X s (t))dt Vip£Ci;(H), 
Jo 
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we infer that R(XI — Lq) is dense in L 2 (H, //) and so the closure L of Lq is cj-m-dissipative 
in L 2 (H, /i) , and it is the infinitesimal generator of the transition semigroup Pt defined by 
(|QD]I . We set N = -L + rl. 

Furthermore, let us assume that 

s(.,zi)€Lip(R ft ), £>^(-,xi)€M(R n ) VxxGFq, (4.25) 

Tr[aa*D 2 X0 g](; Xl ) >0 Vxi G ff , (4.26) 

where (|4.26p is taken in the sense of distributions (or equivalently in the sense of M(W 1 )). 
This implies, as in previous cases, that condition (|2.10j) is satisfied. 

In particular, note that (|4.25p and (|4.26p hold if g = ao9o + a i9ii as i n the example 
mentioned above. Assumptions (I4.25h and (I4.26P imply that 

Ng = -- Tt[aa*D 2 g] - r (x , D xo g) Rn - {xi,D xl g) Hq + rg 
< -r (x ,D X() g) Rn - (xi,D Xl g) Ho +rg, 

hence (Ng) + G L 2 (H,fi), because (x ,D xo g) Rn , (xi,D Xl g) Ho , g G L 2 (H,[i). 

Once again the results established in sections [3] allow us to characterize the price of 
the American option as solution (mild, in general, as the typical payoff function g is not 
smooth) of the backward variational inequality on [0, T] 

^-N<p-Ar g (<p)B0 

with terminal condition (f(T) = g, where M g is the normal cone to 

Kg = {x G H : (f(x) > g{x) fj,-&.e.}. 
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